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ABSTRACT 

We study the dynamical evolution of the young star cluster Arches and its dependence on the 
assumed initial stellar mass function (IMF). We perform many direct iV-body simulations with various 
initial conditions and two different choices of IMFs. One is a standard Kroupa IMF without any mass 
segregation. The other is a radially dependent IMF, as presently observed in the Arches. We find that 
it is unlikely for the Arches to have attained the observed degree of mass segregation at its current age 
starting from a standard non-segregated Kroupa IMF. We also study the possibility of a collisional 
runaway developing in the first ~ 2 — 3 Myr of dynamical evolution. We find that the evolution of this 
cluster is dramatically different depending on the choice of IMF: if a primordially mass segregated 
IMF is chosen, a collisional runaway should always occur between 2 — 3 Myr for a broad range of initial 
concentrations. In contrast, for a standard Kroupa IMF no collisional runaway is predicted. We argue 
that if Arches was created with a mass segregated IMF similar to what is observed today then at the 
current cluster age a very unusual, high-mass star should be created. However, whether a collisional 
runaway leads to the formation of an intermediate-mass black hole (IMBH) depends strongly on the 
mass loss rate via winds from massive stars. Growth of stellar mass through collisions can be quenched 
by strong wind mass loss. In that case, the inter-cluster as well as intra-cluster medium are expected 
to have a significant Helium enrichment which may be observed via Helium recombination lines. The 
excess amount of gas lost in winds may also be observed via X-ray observations as diffused X-ray 
sources. 

Subject headings: methods: N-body simulations, methods: numerical, globular clusters: individual 
(Arches), stellar dynamics 



1. INTRODUCTION 

One of the biggest uncertainties in star cluster evo- 
lution studies lies in determining the true initial stel- 
lar mass function (IMF). Traditionally all studies fo- 
cusing on the dynamical evolution of star clusters as- 
sume an initially fully formed cluster with a given num- 
ber of single and binary stars, all at zero age main 
sequence (ZAMS) at t = and no gas. It is also 
often assumed that the IMF is a standard power-law 
(or power-laws with different indices in different mass 
ranges) wi th no primordial radial variation in t he clus- 
ter ( e.g., ISalpeterl fl959t IMiller fc Scalol H979t IKroupal 
120011 ). However, star formation simulations indicate that 
neither of the above simplifications may be valid. For 
example, star formation takes place over a time-period 
and the formation timescale depends on many different 
physical processes including the Jeans mass of the pro- 
tostellar clouds , radiative feedback and turbulence (e.g., 
iTan et aLll2006f ). Moreover, depending on these detailed 
physical processes stars of different mass ranges may 
form over different t imescales (e.g., iDib et all l2007at 
iKrumholz et al.ll2009f ). In addition, deviations from the 
standard MFs are observed in many cluste rs both at 
the h igh mass and the low mass ends (e.g., Elmegrecn 
120041 ). In particular at the high mass end observed MFs 
are generally top-heavy compared to standard power- 
laws in dense clust er cores like the Arches cluster and 
starburst regions dStolte et all 120021 : iKim et all 120061 : 
iDabringhausen et al.ll2009f ) . Mass segregation is also ob- 



served in old globular clusters (e.g.. iSubramaniam et al 



iHillenbrand fc Hartmannl 119981 : iBaumgardt et al 
Whether the observed mass segregation is im- 
printed from the star formation epoch or it is attained 
via dynamical evolution of the cluster at a later time is 
still an open question. Although a conclusive general an- 
swer to this question does not exist, it has been shown 
in the case of the very young Trapezium cluster in the 
Orion nebula that the observed degree of mass segrega- 
tion cannot be explained simply via dynamical mass seg- 
regation and a preferential formation of high-mass stars 
near the center is necess ary to explain the observations 
(jBonnell fc Daviedfl998f) . 

In order to bypass the above mentioned assumptions 
a detailed model for the star formation processes dur- 
ing the formation of the cluster is needed. Interstel- 
lar clouds contract and fragment due to density per- 
turbations to form pre-stellar cores (PSC) that col- 
lapse on timescales th at depend on t heir Jeans masses 
to form stars (e.g., iNakanol I1966D . Depending on 
the compactness of the assortment of these PSCs they 
can merge with other PSCs and grow before they 
collapse to form stars (e.g., iSilk fc Takahashl 119791 : 
lElmegreen fc Shadmehrl l2003t Ipib et al l 12007b]) . Us- 
ing this simple model IDib et al.l (|2007aD have attempted 
to eliminate the above-mentioned assumptions regard- 
ing cluster IMFs by building an IMF from a distribution 
of PSCs in a molecular cloud. In their simplified model 
they assume that each PSC collapses to form a single star 
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of mass proportional to the parent PSC. Comparing the 
collision timescale between the PSCs with the collapse 
timescale of the PSCs they show that near the central 
regions of a dense cluster the PSCs collide and grow in 
mass before they can collapse to form a star. Thus, near 
the core where the number density of the PSCs is sig- 
nificantly higher than that in the halo, it is more likely 
to form stars that are heavier than stars formed from 
PSCs in the halo. Using the Arches cluster as a particu- 
lar example they show that their semi-analytical model 
naturally results in a mass segregated cluster similar to 
the Arches cluster and the observed top-heaviness near 
the center is thus explained as imprinted from the star 
formation epoch. 

The effects of primordial mass segregation on the 
global dynamical evolution of clusters have been explored 
previously in a parametric way showing interesting dif- 
ferences in the global propertie s of those clusters and 
their dynamical ev olution fe.g.. iBaumgardt et alj [20081 : 
iGiirkan et al]|2004l) . The observed top-heavy MF in the 
central reg ions of the Arches cluster (e.g., iKim et all 
120061 120071 ) thus makes it an interesting candidate for 
such studies. If the observed radially dependent top- 
heaviness in Arches MF is indeed primordial, then this 
is a direct measurement of the degree of mass segrega- 
tion at formation of this cluster, thus it eliminates the 
necessity for parameterization. 

The Arches cluster is also unique for the follow- 
ing reasons. The cluster has a p rojected distance of 
~ 30 pc from the Galactic center flNaeata et al.l [19951 : 
iPortegies Zwart et al.1 [20021: IKim et al J 120061 120071) . Es- 
timates from detailed proper motion observations with a 
baseline of ~ 4yr indicate that the real G alactocentric 
dista nce of the Arches cluster is » 27 pc (jStolte et al.l 
2008). Thus the Arches is the closest massive cluster 
from the Galactic center. The proximity of the cluster 
from the Galactic center necessitates this cluster to be 
compact to be saved in th e face of possible early tidal 
disruption. iSerabvn et alj (|1998f ) a decade ago found 
> 100 O stars within 0.6 pc of the cluster, indicating 
Arches to be the densest young cluster in the Galaxy. 
The high estimated total cluster mass (~ 10 4 — 10 5 M ), 
and the high luminosity 10 s Lq thus make this young 
clust er a power house c luster for Galactic standards 
(e.g. . ISerabvn et aT1ll998f K The central st ellar density of 
the Arches cluster is ~ 10 5 Mfn/p c 3 (e.g.. ISerabvn et all 
1 19981: IPortegies Zwart et al.ll2007h . which is typical only 
for some old Galactic globular clusters. The core radius 
(r c ) of the Arches cluster is < 0.2 p c (e.g., iFiger et all 
[19991: IStolte et al.ll200"i IStoltd 12006). The high central 
density and the high mass thus makes the Arches cluster 
a very good candidate to study for possibility of colli- 
sional runaway. The Arches cluster is also very young, 
the c luster age (t c i) is only 2 ± 1 Myr (e.g.. IStolte et alJ 
2008) . The timescale for a successful collisional runaway 
is only ~ 3 Myr since past that age, the cluster is depleted 
of its high-mass stars due to Supernova explosions (e.g., 
iFreitag et al. 2006a, 2007). Hence, if conducive for a col- 
lisional runaway, the Arches cluster may be undergoing 
that process at the present time. 

Two main scenarios to create the observed top-heavy 
radially dependent MF have been suggested in previous 
studies: 1. The observed MF is primordial, as men- 



tioned earlier (e.g., iDib et aLll2007a|) . 2. The observed 
MF is a result of dynamical evolution starting from a 
standard IMF with no prim ordial radial variation (e.g., 
IPorteg ies Zwa rt et al.1 12007). Using numerous numeri- 
cal simulations we study the dynamical evolution of the 
Arches cluster assuming the observed MF was indeed pri- 
mordial. We compare the results from these simulations 
with the results from another set of simulations where the 
initial cluster properties are kept unchanged but a stan- 
dard IMF with no radial variation is used instead. We 
investigate whether it is possible to attain the observed 
degree of mass segregation starting from an unsegregated 
cluster. We verify that the MFs and cluster parameters 
of the simulated clusters at t c i w 2 Myr are consistent 
with observations. We then investigate the differences 
in the dynamical evolution of the simulated Arches-like 
clusters with the two above-mentioned choices of IMFs 
focusing on the possibilities of a collisional runaway. Fur- 
thermore, we study the final fate and discuss the possible 
observational signatures of a collisional runaway depend- 
ing on the stellar wind prescription. The arrangement of 
this paper is as follows. In 5J5]we describe our simulations. 
In fJ3]we summarize our key results. This is followed by 
2] where we discuss the key observable signatures of a 
collisional runaway. We summarize and conclude in SJSJ 

2. INITIAL CONDITIONS 

We simulate clusters consisting initially of 65536 stars 
with po sitions and velocities chosen according to a King 
profile (|Kind[l966l ) with dimensionless central potential 
Wq in the range from 5 to 9. The initial virial radius 
(r v ) of the cluster is set to 0.5 pc, which results in simu- 
lated clusters that have characteristic radii, e.g., r c and 
tidal radius (r t ) at w 2 Myr close to the observed values 
for the Arches cluster (Table [TJ see also iJT]for observed 
values). Our simulated clusters ar e initially Roche lobe 
filling . Numerical simulations by iGieles fc Baumgardtl 
(2008) show that dynamical evolution and tidal trun- 
cation makes initially underfilled and overfilled clusters 
Roche lobe filling and Roche lobe limited at timescales of 
the order of their half-mass relaxation time (t r j t ). Since 
the typical initial t r ^ for our simulated Arches-like clus- 
ters is much lower (~ 10 3 yr) than the current age of 
Arches, the above assumption is justified. 

All the simulated clusters are assumed to be in a circu- 
lar orbit around the Galactic center. Although the clus- 
ter orbit may be eccentric, the variations of its distance 
to the Galactic center are not ex pected to be significant 
for this study ijStolte et al.ll2008f) . We simulate all clus- 
ters at a fixed Galacto-centric distance since given the 
present day mass and Galactic location of the Arches 
cluster, inspiral timescale from a larger distance is more 
than an order of magnitude longer than its estimated 
age (jGiirkan fc Rasidl2005l: TStolte et al.ll2008f ). One may 
argue however, that Arches ini tially was a lot more mas- 
sive (needs to be > 10 6 M Q , IKim fc Morriil2003D and 
has lost most of its mass during its evolution and in- 
spiral and is now observed at its current position just 
before complete disruption. We find in our simulations 
that the total mass loss from escaping stars through the 
tidal boundary and from stellar evolution typically is 
only ~ 0.04 of the initial cluster mass during the first 
2 Myr. Since, at the current position of Arches the tidal 
effects of the Galaxy is the maximum if inspiral scenario 
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Fig. 1. — Comparison of the Seg IMF within 2r c from the 
cluster center of mass with other standard IMFs and the ob- 
served data points. Black dotted li n e is t he standard Salpeter 
IMF, red dashed line is the IKroupal J20011') IMF, and The Blue 
dots and the error-bars are e xtracted us i ng AD S's Dexter applet 
UDemleitner et al.l [20011) from IDib et al.l I2007af ) showing the ob- 
served MF within 2r c in Arches. The blue solid line is an approxi- 
mation of the observed MF within 2r c in the Arches cluster using 
broken power-laws of the form n oc M~ a , where a changes in dif- 
ferent mass ranges. The a values are given in the form Mi(a)M2 
denoting the value of a in the mass range Mi to M2 in M©: 
0.1(1. 3)0.5(2.3)1.0(2.04)3.0(0.)15(1. 72)100. For r > 2r c , the IMF 
Seg is standard Kroupa 2001. 

is to be believed, this mass loss rate is an upper limit. 
Thus it is unlikely that Arches was initially more massive 
thanl0 6 M Q . 

For each cluster profile we perform two sets of simula- 
tions using different IMFs to assign masses to stars. For 
the set of simulations with init i ally n on-segregated clus- 
ters we use the IMF in IKroupal (|2001[ ) for the whole clus- 
ter (henceforth, Nseg). For the second set we choose a 
radially dependent IMF (henceforth, Seg), with a central 
top-heavy mass function fo r r < 2r e , very similar to the 
observed one of iKim et alj (pOOrl Figure EJ and a stan- 
dard IKroupal (l2001n " IMF outside. Initial cluster proper- 
ties and key results are summarized in Table [TJ In order 
to make sure that our results are robust against statisti- 
cal fluctuations we repeat each simulation 2 times varying 
the random seed for generating the positions and veloc- 
ities for each cluster with the same IMF and initial Wo. 
In order to advance the cl usters in time we use the direct 
TV-body code NBODY4 (|Aarsethl [2003] ) which includes 
the tidal effects of the Galaxy, the ev olution of stars us- 
ing th e integrated SSE/BSE package (jHurlev et alJl2000L 
2002) , and makes use of the GRAPE6 hardware to speed- 
up the computation of the gravitational forces between 
stars. The evolution of each cluster is followed for 5 Myr, 
which is longer than the approximate time when the first 
supernova explosions occur, turning the most massive 
stars into black holes and consequently ending any possi- 
ble run away-growth (|Freitag fc Benz 2002; Freita g et alJ 
l2006bl lal). We do not have primordial binaries in our 
simulated clusters, although binaries created via three- 
body interactions are automatically included. The ini- 
tial ~ 10 Myr evolution of a star cluster is dominated 
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Fig. 2. — MFs of the simulated clusters for runs SW8-1 (top panel) 
and NSW8-1 (bottom panel) at 2.0 Myr. The red (dashed) and the 
black (solid) histograms show MFs outside and inside r = 2r c , 
respectively . The blu e dots show the present day observed MF, 
taken from Dib ct al. (2007a) using the Dexter application from 
ADS. The simulated histograms are normalized so that the area 
enclosed is unit y. The observed p oints have the same normaliza- 
tion as used by IDib et all 12007al . see right labels on the vertical 
axis). Both simulations show some degree of mass segregation. 
The primordially unsegregated cluster NSW8-1 do not attain the 
observed degree of mass-segregation at t c i = 2 Myr. In contrast, 
the simulated MF matches the observed MF reasonably well for 
run SW8-1 with the primordially mass-segregated cluster with the 
Seg IMF. 

by stellar evolution mass loss and two-body relaxation. 
Strong interaction mediated by binaries become a more 
important factor only at later times so far as the global 
evolution of the cluster is concerned. Including binaries, 
however, may decrease the mass segregation timescale 
and increase collision rate to some degree. Nevertheless, 
since the binary fraction in Arches is not known obser- 
vationally, we choose not to consider primordial binaries 
rather than making a wild guess; the additional benefit 
is a significant reduction in computational cost. 

3. RESULTS 

In this section we present the results from our numeri- 
cal simulations. We show in detail the results for clusters 
with Wq = 8 and we summarize more briefly the results 
for all other runs. 

3.1. Evolution of the Mass Function and Observational 
Constraints 

Here we show, how well the initially segregated (Seg) 
and unsegregated (NSeg) cluster models can reproduce 
the observed level of mass segregation at the current age 
(2 Myr) of the Arches cluster. As a condition for consis- 
tency, the simulated clusters should have a present day 
MF consistent with the observed MF at the age of the 
Arches cluster. We compare the MF from the simulated 
clusters with (Seg) and without (NSeg) mass segregation 
at t c i = 2 Myr with the observed MF. Figure [2] shows 
the MF within r = 2r c for runs SW8-1 and NSW8-1 at 
td = 2 Myr. Each mass bin contains 50 stars so that 
the Poisson error is equal in each bin, and the MF is 
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TABLE 1 
List of simulations 



Name 


Al [10 


MqJ 


King 


IMF 




V 


(P c ) 




Mmax (,M© j L 


#Coll 




Myr 


2 Myr 


Wo b 




r c 


Tt 










SW9-1 


ft o 


k n 

o.y 


y 




n i 7 _L n nc;f 
U. 1 ( ± U.UO 


in 7 ±n t 
iy. ( ±0.( 


n 7Q _l n riQ 
U. f o ± U.Uo 


i 1 a J_ n nft 
1. 14 it U.UO 


QQ7 OC^O 


iy 


SW9-2 


ft Q 
. O 


ft i 

0. 1 


y 




n i o _L n n/i 
U.1Z ± U.U4 


1 t I n ft 
1 / .o ± U.u 


n ftft J_ n m 
U.OO ± U.U1 


i 1 1 _L n n n7 
1.1 / it U.U.U / 


Q Kft K 

oOO, - & 


y 


CM.TS3 1 


O 

y.z 


e s 
o.o 


D 
O 




n ii 4_n nQ 
U.ll it U.Uo 


io.z ± u. / 


n fto 4- n m 
U.OZ it U.U1 


n 7/i -U n n/i 

U. ( 4 it U.U4 


QOft QOft 

oUO, ozu 


oU 


SW8-2 


n n 
y.U 


Q ft 


Q 




n n? _l n no 
U.U ( ± U.UZ 


n * i n k 
y.4 ± U.O 


U.4o it U.U1 


n Qft _L n no 
U.oo it U.Uz 


i nn7 ft/i i 
!Uy<, 04 f 


•1 T 7 
4 / 


SW7-1 


n o 

y.z 


o n 

o . y 


T 
I 


Seg 


n i n J_ n no 
U.1U ± U.UZ 


n ft _L n q 
y.D ± U.o 


n /in x n m 

u.4y it u.ui 


n a ft J_ n no 
U.40 it U.Uz 


ftCQ 07Q 


o o 


SW7-2 


inn 


("1 T 

y. i 


7- 
( 


n i q J_ n n^i 
U. lo ± U.U4 


1 n q _l n q 
lU.o ± U.o 


n c^o J_ n m 
U.oZ it U.U1 


n ^ tn J- n no 
U.OO it U.Uz 


7HQ /lOO 

(Uo, 4zo 
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SW6-1 


IOC 

IZ.o 


1 O A 

1Z.4 


c 
O 




n i n J_ n no 
U. 1U ± U.UZ 


Q A _i_ n o 
0.4 ± U.Z 


n yi ft _L n m 
U.4D it U.U1 


r\ An J- r\ no 
U.4U it U.Uz 


i /i n 1 1 7 
14U, 11/ 


lo 


SW6-2 


ion 

iz.y 


1 O e; 
1Z.0 


a 





n 1 1 _i_ n m 
U.l * ± U.U1 


q ? I n o 
o. f ± U.Z 


n a i J_ n m 
U.4 ( it U.UI 


r\ a a J- r\ ni 
U.44 it U.UI 


QQO 1/1/1 

ooz, 144 


Oft 


SW5-1 


ICQ 


10. o 


- 

O 




n 1 n _i_ n m 
U. ly it U.U1 


? ft i n o 
I .0 ± U.Z 


n a q J- n m 
U.4o it U.UI 


n /I n _L n ni 
U.4U it U.UI 


i oft 1 oc; 
1Z0, 1Z0 


"1 1 
11 


SW5-2 


15.5 


15.0 


5 




0.20 ±0.01 


7.4 ±0.2 


0.43 ±0.01 


0.39 ±0.01 


156, 107 


13 


NSW9-1 


3.23 


3.18 


9 




0.04 ±0.01 


12.0 ±0.3 


0.55 ±0.01 


0.8 ±0.2 


339, - 


14 


NSW9-2 


3.20 


3.18 


9 




0.04 ± 0.04 


11.9 ±0.2 


0.55 ±0.01 


0.2 ±0.2 


238, 104 


15 


NSW8-1 


4.2 


4.1 


8 




0.18 ±0.09 


12.2 ±0.8 


0.55 ±0.02 


0.57 ±0.06 


136, 131 


14 


NSW8-2 


3.2 


3.2 


8 




0.06 ± 0.02 


8.9 ±0.2 


0.47 ±0.01 


0.28 ± 0.06 


135, - 


8 


NSW7-1 


3.15 


3.14 


7 


NSeg 


0.08 ±0.04 


7.6 ±0.1 


0.4 ±0.1 


0.4 ±0.1 




8 


NSW7-2 


3.23 


3.21 


7 


0.11 ±0.03 


7.6 ±0.3 


0.436 ±0.008 


0.24 ± 0.05 


130, - 


5 


NSW6-1 


3.3 


3.2 


6 




0.08 ±0.02 


7.3 ±0.1 


0.426 ±0.004 


0.26 ± 0.02 


116, - 


2 


NSW6-2 


3.2 


3.2 


6 




0.11 ±0.02 


6.8 ±0.1 


0.412 ± 0.003 


0.37 ± 0.04 




3 


NSW5-1 


3.2 


3.1 


5 




0.15 ±0.05 


15.3 ±0.5 


0.62 ±0.01 


0.93 ± 0.07 


233^ 112 


16 


NSW5-2 


3.2 


3.2 


5 




0.08 ±0.03 


6.9 ±0.1 


0.42 ±0.01 


0.12 ±0.02 


147, 104 
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a Bound mass at the age (t c i) of the simulated cluster. 

b King profile HKing||1969) 

c The range in mass is [0.1 - 100] M 

d at t cl = 2 Myr 

c Highest mass reached for the two most massive stars in the simulation. 

f All quoted errors are ltr error-bars from statistical fluctuations between the snapshots over which the cluster properties are averaged 
for the the simulated clusters. 

g Mass does not grow above 100 Mq within 4 Myr 

normalized such that its integral is unity. As it can be 
seen, while the evolved MF in the inner region of the Seg 
cluster is in good agreement with observations, this is 
not the case for the NSeg cluster. For the NSeg cluster 
the slope at the high mass end increases only marginally 
(from -2.3 to -1.6). While for the Seg cluster the MF 
becomes flatter, but remains still compatible with obser- 
vations. Same results are obtained when the MFs are 
compared at t c i w 3 Myr. This indicates that the de- 
gree of mass segregation observed in the Arches cluster 
is not likely to result simply from the cluster's dynamical 
evolution during its lifetime without any primordial mass 
segregation. It also indicates that, the MF at t c i = 2 Myr 
for the simulated cluster is consistent if the observed MF 
of the Arches cluster truly is i mprinted fr om its star for- 
mation epoch as sugg ested by|Dib (2001). 

We should note that for the NSeg IMF, a much higher 
degree of mass segregation can be achieved if stellar evo- 
lution is not included in the simulations. We find that if 
we turn off stellar evolution, the degree of flattening in 
the MF of the simulated cluster at ~ 2 Myr increases sig- 
nificantly compared to that with stellar evolution (Figure 
[3]) . Although the wind mass loss from the high-mass stars 
due to stellar evolution is negligible (< 4%) compared to 
the total cluster mass in these simulations during the first 
2 Myr, this mass is lost from the most massive stars in 
the cluster, which in turn reside near the center due to 
mass segregation. Thus this mass is lost from the deepest 
part of the cluster potential and can have a perceptible 
effect on the global properties of the cluster. From our 
results we conclude that the mass loss via stellar evolu- 
tion from the massive stars can inhibit mass segregation 
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Fig. 3. — MFs of the simulated clusters for runs NSW8-1 with (top 
panel) and without (bottom panel) stellar evolution at 2.0 Myr. 
The red (dashed) and the black (solid) histograms show MFs out- 
side and inside r = 2r c , respectively. The blue dots show the 
present day observed MF, taken from (Dib ct al. 2007a) using the 
Dexter application from ADS. The simulated histograms are nor- 
malized so that the area enclosed is u nity. The observe d points 
have the same normalization as used by Dib et al. (2007a, right la- 
bels on the y axis). In both cases the initial clusters are identical. 
The simulation including stellar evolution shows significantly lower 
degree of mass segregation than that without stellar evolution. 
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Fig. 4. — The growth of mass in time for the two most mas- 
sive stars in the cluster are plotted for runs SW8-1 (solid lines) and 
NSW8-1 (dashed lines). Each dot denotes the time when the partic- 
ular star changed mass through collision. The area of each dot is 
proportional to the stellar mass. For the Seg IMF a collisional run- 
away is observed starting around 2 Myr. For the NSeg IMF using 
the same initial conditions there is no runaway. 

to a certain extent and simulations done without tak- 
ing that physical effect into account will overestimate 
the amount of mass segregation. Since this effect is only 
predominant for the most massive stars at a given time 
in the cluster, the quantitative effect will depend on the 
number of high-mass stars present in the cluster. De- 
spite a lower total nu mber of initial stars (N ~ 12K and 
simulations in iPortegies Zwart et al.l (|2007f ) had 
more massive stars than our NSeg clusters because of the 
truncated IMF they used. Thus, if stellar evolution is 
included, mass segregation should have been even more 
strongly suppressed in their simulations compared to our 
simulations with the NSeg clusters (Figure [3|) . This fur- 
ther supports the notion that the level of observed mass 
segregation in the Arches cluster cannot have a fully dy- 
namical origin. 

3.2. Choice of IMF and Implications for Collisional 
Runaway 

In general we find that for the clusters with primordial 
mass-segregation (set Seg) the number of collisions is al- 
ways significantly larger compared to the clusters with 
a standard Kroupa IMF without any primordial mass 
segregation (set NSeg; see Table [1]), for example, in the 
case Wq = 8 by more than a factor of two (30 for run 
SW8-1 and 14 for run NSW8-1) within the first 4 Myr. 
This is mainly due to the top-heavyness of the Arches 
IMF within r < 2r c (fJ3J) and the larger overall mass of 
our clusters with primordial mass segregation, resulting 
in a larger number of massive stars with large collision 
cross sections. Consequently, the runaway growth of a 
massive star is much more likely to take place in clusters 
from set Seg. 

As an example, Figure 2] shows the mass of the two 
most massive stars over time for the cluster with the Seg 
IMF (run SW8-1) and the NSeg IMF (run NSW8-1). In 
the Seg case there is a clear evidence of a collisional run- 



away producing two very massive stars. The most mas- 
sive star in this case attains a mass of 800 M Q through 12 
successive collisions. The second most massive star also 
grows significantly and reaches 300 M through 4 suc- 
cessive collisions. A clear sign of a collisional runaway 
is a steep increase of the stellar mass through successive 
collisions. Figure 2] clearly shows that at i c ; ~ 2 Myr the 
most massive star enters a runaway regime. The whole 
process of the runaway takes place between ~ 2-4 Myr. 
In case of NSeg IMF the two most massive stars expe- 
rience only two to three collisions and there is no run- 
away happening. We obtain qualitatively similar results 
when changing the random seed to generate other equiv- 
alent initial conditions. For instance, in the run SW8-2 
with Seg IMF we see a runaway, where the most massive 
star grows to ~ 10 3 M© within the first 4 Myr starting 
from about 90 M , through 17 successive collisions. The 
corresponding run with the NSeg IMF, NSW8-2, did not 
exhibit any runaway growth (Table [1]) . A similar result 
is obtained for simulated clusters with other Wq values. 
For example, Figure [5] shows the growth of the two most 
massive stars in the simulated clusters for runs SW7-1 
and NSW7-1, both clusters with Wq = 7, but the first 
from the Seg and the latter from the NSeg set. Again, the 
run from set Seg shows clear evidence of a collisional run- 
away, whereas the equivalent cluster from set NSeg shows 
none. (For a full list of simulation results with other 
initial Wo values, see Table [T]) Note that, in the colli- 
sional runaway of the stars in our simulations, we find 
that the stars can grow even after 3 Myr via collisions. 
This is due to the collisional rejuvenation prescription 
in BS E, which assumes complete mixing ([Hurley et al.l 
2000). Thus the collisional stars live longer compared 
to their normal lifetime. This assumption is of course a 
simplification and the amount of mixing and thus reju- 
venation depen ds on the details of the kinematics of the 
collisions (e.g., I Sills et al.lfl997l : IGlebbeek fc Polsl[200l 
. However, this assumption gives an upper limit of the 
degree of rejuvenation and lets the stars grow for a longer 
time than with a more realistic rejuvenation prescription. 
In this study for the considerations of a collisional run- 
away we take collisional sequences only upto t c i = 4 Myr. 

A very interesting aspect of these results with the Seg 
IMF is that, many of these runs lead to a double colli- 
sional runaway, one being stronger than the other (e.g., 
runs SW8-1.2, SW7-1.2; see Table [1] for a full list) cre- 
ating two very massive stars (VMS). Double runaways 
have been studied and observed in s ome previous simula - 
tions with primordial binaries (e.g., iGiirkan et al.| [2006). 
however, a double runaway even without primordial bi- 
naries has not been observed before. Due to the flat 
mass spectrum in the high mass end of the Seg IMF, 
there is a larger reservoir of high-mass stars increasing 
the chance of a double runaway. A double runaway does 
not necessaril y mean black hole binaries, as discussed in 
IGiirkan et all (f2006h . since the two VMSs may still col- 
lide with each other at a later time but before compact 
object formation (e.g., SW8-2, SW7-1). The signature of 
a double runaway and production of two VMSs seems 
to be common in all of our simulations with the Seg 
IMF, where the initial concentration is conducive to a 
collisional runaway in particular for runs with King con- 
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Fig. 5. — Same as Figure [4] but for runs SW7-1 and NSW7-1. No 
star in run NSW7-1 grows to > 100 Mq within the first 4 Myr, hence 
are not not shown on this plot. 
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Fig. 7. — Effect of wind mass loss for the run SW8-2. Top panel: 
Solid line shows the mass growth via collision with the standard 
stellar evolution prescription in Nbody4. The circles show the 
time when the collisions happened. The dashed line shows the 
same collision chain, only with the high wind mass loss prescrip- 
tion (Glcbbeck et al. 2009). Following each collision, the collision 
product is evolved and wind mass loss is taken into account. Bot- 
tom panel: The wind mass loss rate as a function of time. 
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Fig. 6. — Masses of the two most massive stars are shown for each 
simulated cluster with the Seg (top panel) and NSeg (bottom panel) 
IMFs as a function of the initial King concentration parameter 
Wq. The filled point and the open point show the highest and the 
second highest mass attained during the first 4 Myr of evolution. 
Cases where the mass is not above 100 Mq are not shown since the 
initial high mass cut-off of the IMF is 100 Mq in our simulations. 
Multiple realizations of the same cluster in each case is simulated 
to ascertain robustness against statistical fluctuations. Note that 
for the Seg IMF with Wo 7 and 8, there is clear indication of double 
runaways. 



centration parameters Wq = 7,8,9 (Figure [6]). 

Similar results are obtained for clusters with other ini- 
tial concentrations or Wq values. The main difference 
is in the final mass of the most massive stars, shown in 
Figure [Sj The final stellar mass is also subject to sta- 
tistical fluctuations to a certain extent. The collisional 
runaway with the Seg IMF and Wq — 6 is weak, whereas, 
for Wq = 5 there is no runaway (Table [J) . For each of 



these runs with the Seg IMF, the corresponding simu- 
lations with the same initial Wq but the NSeg IMF do 
not exhibit any runaway and do not produce the VMSs. 
However, for Wq = 9, both the Seg and the NSeg IMF 
simulations show runaways and the masses of the colli- 
sionally created most massive stars are comparable. In 
the Seg case Wq — 9 shows a weaker collisional runaway 
in our simulated examples. This is because from a King 
model with Wo = 8 to one with Wo = 9, the core radius 
decreases considerably, thus the top-heavy IMF is within 
a smaller volume. As a result the total initial mass as 
well as the number of high-mass stars decrease. From 
our simulations, we find that with the Seg IMF, Wq = 7 
and 8 are optimal for strong collisional runaways. On one 
hand a lower Wo decreases the probability of a runaway 
because of the lower central density. On the other hand, 
a higher W gives the same result although the central 
density is higher, because of the lower overall mass of the 
cluster as well as the lower number of high-mass stars. 

3.3. Enhanced Mass Loss via Stellar Winds 

High-mass stars can lose a significant f raction of 
their mass through s tellar winds (e.g.. iCassinellil 
19791 iChiosi fc Maederl [l98l iKudritzki fc Puis! 12000). 



Glebbeek et all (|2009l ) show that a VMS generated via 



stellar collisions may lose mass through stellar winds 
at such a high rate that its collisional growth may be 
quenched. We explore this possibility, following their 
work, for our simulated clusters SW8-1 and SW8-2, where 
collisional runaways have been observed. Following 
iGlebbeek et alJ (|2009f ) each collision is treated using the 
"Make Me a Massive Star" (MMAMS) soft ware package 
(jGaburov et al.ll2008t lLombardi et al.ll2002D . which gives 
the properties of the collision product. This output is 
then imported into the Eggleton stel lar evolution code 
(|Eggletonlll97illl972tlPols et al.lll995l ) using the method 
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Fig. 8. — Time evolution of the surface abundance of the most 
massive star in the run SW8-2. Surface abundance of Hydrogen 
(H), Helium (He), Carbon (C), Nitrogen (N), Oxygen (0), Neon 
(Ne), and Magnesium (Mg) are shown. Note the very high surface 
abundance of Helium. 

of iGlebbeek et all (|2008t ). This method is repeated for 
every collision i n the collision chain ( for a more detailed 
description see IGlebbeek et al.l 120091 especially Section 
2.1). 

Figure [7] shows the time evolution of the most massive 
star for the run SW8-2 as an example. The solid line in 
the top panel shows the collisional growth of the stellar 
mass in the runaway observed in our run SW8-2, using the 
standard wind prescriptions from BSE/SSE, included in 
NBODY4. The dashed line in the top panel tracks the 
mass of the star a ssuming the wind mas s loss prescrip- 
tions suggested by IGlebbeek et all |2009). The bottom 
panel shows the mass loss rate via stellar w i nds as a func- 
tion of time according to IGlebbeek et all (|2009f ). With 
this wind mass loss prescription, while collisions increase 
the stellar mass, the increased stellar mass in turn leads 
to enhanced wind mass loss. The growth of the stellar 
mass through the chain of collisions is quenched and in- 
stead of the ~ 10 3 M VMS a star of only ~ 30 M Q is pro- 
duced with this wind prescription. At the final stage the 
star is burning Helium in its core and has a high surface 
Helium abundance (X He ,s ^ 0.8) (Figure [5]). We should 
mention that the modeling of stellar winds from high- 
mass stars is quite uncertain, and various treatments of- 
ten come from extrapolations from observed stars at a 
lower mass. Depending on the prescription wind mass 
loss rates may vary by orders of magnitude. If indeed 
the high wind mass loss prescription is correct, then a 
VMS will not be produced even when a collisional run- 
away takes place. Interestingly, the mass lost via stellar 
winds due to the collisional growth of the stars will be ex- 
tremely Helium rich which may provide observable signa- 
tures (more on this in H4.3[) . Similar results are obtained 
from the collision chain in run SW8-1. 

4. POSSIBLE OBSERVATIONAL SIGNATURES OF A 
COLLISIONAL RUNAWAY 

Here we discuss the possible observable signatures of a 
collisional runaway in a cluster. In the previous section 



Fig. 9. — The evolution of the most massive star created in run 
SW8-2 in the luminosity (L) vs effective temperature (T e ff) plane. 
Solid red line shows the evolution of the most massive star for the 
enhanced wind mass loss prescripti on. The dotted line (blue ) show s 
the Humphreys-Davidson limit (Humphreys & Davidson 1979). 
The observed luminous stars are also plotted for reference. The 
filled black squares denote LBVs, open squares denote observed 
OB s tars in the M agellanic clouds, filled triangles denote stars in 
GCs {Figcr et al. 19981, and the references therein). Two points 
for the Pistol star denote high (H) and low (L) luminosity models 
IjFiger e t al. 1998 ). The blue pen tagons show the observed bright 
stars in Arches (Clark et al. 2009, and the references therein). 

we have shown that depending on the stellar wind pre- 
scription a collisional runaway may create a VMS or the 
mass growth can be quenched through high wind mass 
loss. Here we consider both cases and discuss possible 
observable signatures of both possibilities. 

4.1. Possible Creation of another Pistol Star 

The Quintuplet cluster is very similar to the Arches 
cluster. It is also very near t he Galactic center 
(projected distance < 40 p c; e.g., iNagata et al.l 119901 : 
iPortegies Zwart et alj f2002). The Arches and the Quin- 
tuplet are also simi lar in their metallicities, and mass 
(|Naiarro et al.ll2009l) . The Arches is younger (2 ± 1 Myr) 
than the Quintuplet cluster (~ 4 Myr). The Pistol star 
is one of the most interesting and well studied objects 
in the Quintuplet. T his is one of the brightest stars 
ever observed (Figure l9l lMoneti et al.lll994H C otera et all 
\199& lLang et al.H l997; Fig er et al.lll998T T The inferred 
initial mass of the Pistol star is 200 - 250 M . The Pis- 
tol star is classi fied as a l uminous blue vari able (LBV 
iFiger et al.lll998f ) or a B[e] (|Morris et al.lll996f ) star. The 
stellar surface is extremely Helium rich and it has lost a 
significant fraction of its mass via Helium rich massive 
stellar winds (e.g. Ifiger et all 119951; iCotera et al.l 119961: 
i Timmermann et al.l 119961: lLang et alj 119971 : IFiger et all 
1998). From our results we find that the massive stars 
generated via collisional runaways may be a potential 
channel to produce a massive and Helium rich Pistol 
like star (Figure [8]). Figure [9] shows the evolution 
of the most massive star found in run SW8-2 for the en- 
hanced wind mass loss case as an example. The massive 
star produced through the collisional runaway spends 
a significant period of time (e.g., 1.4 — 2.3 Myr in case 
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of the run SW8-2) above the H umphreys-Davidson limit 
([Humphreys fc Davidson1ll9790 like the Pistol star. Thus 
if the Arches was indeed created with a top-heavy IMF 
at the center as is observed today, then a Pistol-like star 
could be created at some point. However, due to the un- 
certainty in the age estimate of the Arches cluster and the 
onset of the collisional runaway in simulations depend- 
ing on the initial conditions it is not possible to predict 
whether a similar LBV star should already have been 
created in the Arches cluster and could be observed to- 
day. 



4.2. IMBH Progenitors 



The presence 
center of most 



of supermassive black holes at the 
alaxies is now well established (e.g. , 
Richstone et all 119981: iKormendv fc Gebhardtl 120011 : 
Ghez et al.l 120051 ). Extrapolating the relation between 
the host mass and the black hole mass it has been 
postulated that massive clusters may also host black 
holes in the mass range 10 2 - 10 4 M^ (e.g., see re- 
views by Ivan der Mareli 120041 : iRasio et alJ 120041 ) . These 
black holes are much less massive than the supermas- 
sive black holes at the centers of most galaxies but 
more massive than the stellar-mass black holes ex- 
pected as remnants from high-mass stars. Although 
no conclusive observational evidence for an IMBH has 
been obtai ned yet, there are s everal interesting candi- 
dates (e.g.. [Farrell et al.l [20091 : llrwin et alJ [20091 ). The 
main proposed channel of formation for these IMBHs 
is through collisional runaways in dense clusters (e.g. 
Freitag fc Benz||2002t IPortegies Zwart fc McMillanl |2002 
Giirkan et alJ 12004 120061 : iFreitag et all l2006bl lar 
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20063S). It is thus exciting to find young star clusters 
where a collisional runaway may have taken place. 

The simulated clusters with the Seg IMF, in partic- 
ular with Wq = 7, 8, and 9 create VMSs via colli- 
sional runaway. The counterparts of the same clus- 
ters with the more standard NSeg IMF do not produce 
these VMSs. If the winds do not quench the mass 
growth via collisions, the n the VMSs can potentially 
be IMBH progenitors (e.g. lFrverlll999l: [Fryer fc Kalogera] 
2001 1 : IPortegies Zwart fc McMillan! [20021 : iHeger et all 
200a iBelczvnski et al.ll2009h . If an IMBH is indeed pro- 



duced, there can be many observational signatures. Even 
without primordial binaries, at the high densities at- 
tained at the core binaries can form dynamically from 
three-body interactions. IMBHs being a lot more mas- 
sive than the typical stars in its vicinity inevitably ex- 
change into a b inary or a higher order bound system 
([Gill et al.l 120081 ). A stable mass transfer to the IMBH 
from its binary companion can be visible as ultra lu- 
minous X-ray sources. Inspiral and capture of compact 
objects such as stellar mass black holes and neutron stars 
into the IMBH can be a strong gravity wave source de- 
tectable with the Laser Interferometric Space Antennae 
(LISA). The simulated clusters with the Seg IMF, often 
show signatures of two separate runaways, one weaker 
than the other. These runaways may produce two VMSs, 
both of which may create massive black holes. If the 
VMSs do not collide prior to creating the massive black 
holes a massive binary black hole system can be formed. 
If produced, an IMBH-IMBH binar y is a good source fo r 
LISA during the inspiral phase (see lGiirkan et alJ l2006). 



4.3. Enhanced Wind Mass Loss: Helium Enrichment 

If the wind mass loss of the most massive star domi- 
nates over its collisional growth, then the rapid increase 
in mass is inhibited. (Figure [7]). In this case a VMS is 
not produced as a result of the runaway. Instead, a large 
amount of processed stellar material is injected into the 
surrounding interstellar medium, which is, consequently, 
expected to become strongly enriched in Helium. For 
example, the total mass lost via winds from the most 
massive star formed in the runaway in the run SW8-2 is 
~ 1O 3 M , while its Xiie,S can exceed 0.9, resulting in 
ss 500 M Q of Helium ejected via winds. In contrast, if 
the stars, that merged to form the VMS, are evolved in 
isolation, the total mass of ejected Helium would be sig- 
nificantly lower, « 100 M Q , within the first 3Myr. An 
anomalously high Helium enrichment of the surrounding 
interstellar medium may, therefore, provide an observa- 
tional signature for a collisional runaway. 

However, such a signature might not necessarily be 
found inside Arches but rather outside of it, as the ma- 
terial lost through winds has rather large velocities. As- 
suming the terminal flow velocity of the wind from O/B 
stars for solar metallicities is 0.5 times the esca pe speed 
from the surface of the star (jVink et al.ll200l[ ). we find 
that the wind velocity for our collisionally formed VMSs 
is rj 500km/s. As these velocities are much larger than 
the escape speed of even much more massive star clusters, 
it is likely that the material escapes the Arches cluster 
entirely. On the other hand, this material can also slow 
down through collisions with winds from other massive 
stars. Such multiple stellar wind interactions result in 
a complex network of shock compressed, extremely hot 
(10 7 — 10 s K) gas that should then leave the cluster as a 
cluster wind (see, e.g.. lCanto et al.ll2000f). The structur e 
of this wind consists of four zones (Sili ch et al.ll200l : 
a star cluster region filled with a hot X-ray plasma, 
an X-ray halo with decreasing temperature profile, a 
line-cooling, recombining zone, and a region of recom- 
bined gas, exposed to the UV radiation from the cen- 
tral star cluster. Indeed diffuse X-ray bright sources 
have been observed within Arc hes as well as extending 
beyond the cluster rad ius (e.g. lYusef-Zadeh et al.|[2002t 
Rockefeller et al. 2005|). This extremely hot gas may also 
be confined due to ram pressure exerted by a molecular 
cloud or other external medium surrounding the cluster 
([Rockefeller et al.l[27)05h . 

Helium abundances are usually determined 
through measuring Helium recombination lines (e.g., 
iPeimbert et "ail [2000 ). Thus a strong Helium enrichment 
can then only be detected outside the X-ray halo region, 
in the line coo li ng zo ne. Based on detailed calculations 
of ISilich et all (pOOl . this region is outside a radius 
of 3-7 pc from the Arches cluster. Other signatures of 
a collisional runaway might be provided by the wind 
structure itself. This is because the extension and 
physical parameters of the diffuse X-ray emitting region 
are mainly dependent on the sum of the mass-loss 
rates of the individual sources. Given that for the case 
of a runaway this mass loss rate can be significantly 
increased, it may strongly influence the appearance of 
the cluster winds. However, such imprints might be 
very diffucult to detect. This is in part because our 
rates for the most massive star are highly time- variable, 
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varying by two orders of magnitudes within a few 
10 5 yr, while particularly large values in excess of the 
total mass injection rate for all other Arche s stars of 
7 w xlO^Moyr" 1 (|Stevens fc Hartwellll2003l ) are only 
attained during an even shorter timescale. Another 
reason is that the Arches cluster wind may be confined 
by cooler surrounding i nter stellar material, a n effect 
that is hard to quantify (|Rockcfcllc r et al.ll2005ft . 

Instead of searching for imprints of a runaway in the 
global cluster wind structure, signatures may be found in 
the properties of individual objects, as stellar radii and 
m ass loss rates of runaw ay stars can be unusual (Figure 
IGlebbeek et al.l l2009) . One way to directly detect and 
constrain the mass loss rate, the effective radius of the 
star, and the wind velocities is by exploiting the s o called 
Baldwin effect (|Baldwidll9"77tlMorris et al.lll993h . which 
is an anticorrelation between the line luminosity of one 
strong emission line, such as Hell (ofte n found in WR 
winds ), and the continuum luminosity. Ivan Gent et al.l 
(2001) show that the ratio of line and continuum lumi- 
nosity the equivalent width W, depends on the mass loss 
rate, M, the terminal velocity of the wind, Voo, and the 
effective radius , i?*, as W oc M 2 v^R~ 3 . When the 
distance of the star is known, and the temperature can 
be estimated, observation of W can constrain the wind 
mass loss rate, wind velocities and the effective radii of 
the stars. 

5. SUMMARY AND CONCLUSION 

We carried out many iV-body simulations with initial 
conditions chosen su ch as to resem ble the Arches clus- 
ter. Using a standard iKroupal (|2001f ) IMF with no radial 
variation (NSeg) and a radially dependent IMF that cor- 
responds to the level of mass segregat ion inferred from 
observations (Seg; e.g. lDib et al.ll2007al ). we study possi- 
ble consequences and discuss potential observable signa- 
tures to distinguish between the two cases. 

We show that the degree of mass segregation via dy- 
namical evolution of the simulated clusters during the 
first 2 Myr is minor both for the Seg and the NSeg IMF 
(Figure^]). Thus, while runs with the Seg IMF repro- 
duce well all parts of the currently observed slopes in the 
MF at w 2 Myr, those with the NSeg IMF do not pro- 
duce enough mass segregation within the Arches lifetime 
to match the observed MF as well (Figure [2]). There- 
fore, our results indicate that the present day observed 
degree of mass segregation in the Arches cluster cannot 
be a result of dynamics only and at least some degree of 
primordial mass segregation is required. 

We further show that the choice of IMF changes the 
overall dynamical evolution significantly for Arches like 
clusters over a range of initial Wo. For the Seg IMF 
there is clear evidence of collisional runaway, producing 
stars up to ~ 10 3 M© over a range of initial concentra- 
tions (Table [1] Figures El EJ)- In contrast, with the 
more standard NSeg IMF no collisional runaway is ob- 
served with the same initial conditions. We find that 
the sequence of collisions in the simulations with suc- 
cessful collisional runaway starts at around 2 Myr. Note 
that t he current estimate d age of the Arches is 2 ± 1 Myr 
(e.g., iStolte et al.ll2002[) . Hence, if the observed top- 
heavy MF within r ~ 2r c is indeed primordial, then it 
is possible that a collisional runaway has started or is 
bound to happen. Unfortunately, the uncertainty in the 



age estimation prevents us from predicting whether a col- 
lisional runaway should already have taken place in the 
Arches. 

An interesting aspect of the collisional runaways ob- 
served in this study is the general tendency of a dou- 
ble runaway in cases where the collisional runaway is 
strong (e.g., runs SW8-1.2 and runs SW7-1.2, see Table 
[T]). Double collisional runaways have been observed and 
studi ed in numerical simulations before (|Gurkan et al.l 
2006). However, those simulations concluded that pri- 
mordial binaries are a necessity for a double collisional 
runaway. It is very interesting that with the primordially 
mass-segregated clusters even without any primordial bi- 
naries double collisional runaways can take place. Due 
to the flatness in the high mass end of the IMF in this 
case, there is a larger reservoir of high-mass stars. Thus 
the probability to start a second runaway is higher. The 
ultimate fate of the smaller runaway is dependent on the 
statistical fluctuations of the simulation (see . How- 
ever, it is possible that the VMSs do not collide prior to 
compact object formation and both stay bound to the 
cluster after their respective supernova explosions (e.g., 
run SW8-1). In such a case they may form IMBH bina- 
ries or coalesce at the center of the cluster. If possible to 
form, this can create a very strong gravity wave signal 
for LISA (|Gurkan et al.ll2006h . 

If the mass loss from stellar evolution driven winds is 
too high then the col lisional mass growth is quenched 
(jGlebbeek et al.ll2009f ). By evolving the collision prod- 
ucts separately following each successive collision in a 
sequen ce of runaway collision s according to the prescrip- 
tion in lGlebbeek et alj (|2009T) we find that the VMSs are 
not produced as a result of a collisional runaway. In- 
stead, a Helium star of a few tens of a Solar mass is 
created (Figure [7]). The mass lost via these winds is ex- 
tremely Helium rich, with X}je,s reaching beyond 0.9. In 
this case the cluster as well as its surrounding medium 
will be enriched by this Helium rich gas, which could 
be observed through recombination at distances beyond 
3 — 7pc from the Arches cluster. 

Further indications of a collisional runaway might be 
provided by the runaway object itself, as such an ob- 
ject has a rather unusual radius and mass loss rate com- 
pared to ordinary massive stars. Constraints on the 
radius, mass loss rate, and the wind speeds may be 
obtained using the B aldwin effect (e.g.. iBaldwinl fl977l : 
Ivan Gent et alJfeOOlf ). Furthermore, during the evolu- 
tion of the runaway star, a Pistol like stellar object may 
also be produced ( fc|4.ip . However, since the duration 
for which the star attains Pistol like high luminosities is 
short, it is not possible to predict whether a Pistol like 
star should already have been created in Arches given 
the uncertainties in the age estimation and the statisti- 
cal fluctuations in the exact onset and the details of the 
simulated collisional runaways. 

Nevertheless, we should point out that the above recipe 
to follow the runaway star for these high stellar winds is 
not completely self-consistent. Each collision in a clus- 
ter is dependent on the masses and radii of the collision 
progenitors. Since, severe wind mass loss from the stars 
changes the masses of the stars significantly, the chain 
of collisions in the runaway may also be different. How- 
ever, this model still provides us with an estimate of the 
total mass loss from the VMS via stellar winds, its stel- 
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lar properties, as well as the composition of the ejected 
material if a runaway did indeed happen. 
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